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ABSTRACT 

Various theoretical arguments make exotic heavy neutral weakly interacting 
fermions, particularly those predicted by supersymmetry theory, attractive 
candidates for making up the large amount of unseen gravitating mass in galactic 
halos. Such particles can annihilate with each other, producing secondary particles 
of cosmic-ray energies, among which are antiprotons, positrons, neutrinos and y- 
rays. Spectra and fluxes of these annihilation products can be calculated, partly 
by making use of e + e" collider data and QCD models of particle production derived 
therefrom. These spectra may provide detectable signatures of exotic particle 
remnants of the big bang. 


INTRODUCTION 


Astronomers have learned that there is much more to the universe than meets the 
eye*. The discovery that the rotation curves of spiral galaxies are quite flat to 
radii as far out as can be observed indicates the presence of dark matter comprising 
most of the mass in galactic halos and which could also dominate the mass of the 
universe 2 . This dark matter is most likely of non-baryonic form . It may be made 
up of exotic new particles, perhaps such as those predicted by supersymmetry 
theory. The large body of literature on various aspects of the dark matter problem 
and dark matter detection 5 have recently been reviewed. . . 

Possible candidates for dark matter are exotic, stable, weakly interacting 
particles, which we hereafter call x particles. Among these, the neutralinos 
predicted by supersymmetry theory should be left over as relics of the early big 
bang in cosmological ly significant densities ’ . Such particles can annihilate o 
produce baryons of cosmic-ray energies as well as other familiar particles w ose 
decays will produce cosmic-ray leptons (electrons, positrons and neutrinos) and also 
cosmic y-rays. The fluxes and energy spectra of these ''ordinary particles 
resulting from the annihilation of various x particles can be calculated. Ihe 
annihilation process will produce a characteristic high energy cutoff in these 
spectra at the rest mass energy of the x particle, typically -5 to -30 GeV. e 
Lund Monte Carlo 5 simulation technique used by particle physicists in comparing wijh 
e + e _ collider data provides a powerful tool for predicting xx annihilation spectra 
* 3 , since the process of quark-antiquark pair production is involved in both 
cases. We will present and discuss detailed energy spectra from the annihilation of 
x particles in the Galaxy and show how the features of these spectra may help lead 

to indirect evidence of exotic dark matter particles. 

The liqhtest of the neutral supersymmetric particles, designated the LSP, would 
be stable by virtue of a new natural conservation law called R-parity conservation, 
since it is the lightest state with odd R-parity. The mass eigenstate neutralino 
can be a pure state, but is more generally a superposition of the higgsino (h), 
"photino" (v) and "zino", which are mixed by gauge and supersymmetry breaking. 
There is also the possibility that the x particle could be a heavy Dirac neutrino. 
However, unlike the case with the LSP, there is no natural way of forbidding the 
decay of a "conventional" heavy neutrino. There is also experimental evidence 
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against dark matter heavy Dirac neutrinos.^ 

DARK MATTER ANNIHILATION 

The basic physical processes involved in the xx annihilation process are 
portrayed schematically in Fig 1. In the first stage of the annihilation (Stage "a" 
in Fig. 1), x particles annihilate via a propagator a ( e.g. , a scalar fermion in the 
case of photinos and a Z particle in the case of higgsinos). These supersymmetric 
neutralinos are Majorana fermions, i .e. , they are self annihilating. Then, (b) the 
annihilation usually produces a quark antiquark pair ( xx -*■ qq ). This is followed 
by (c) a cascade of the quarks down to the mass shell by successive gluon emission, 
producing a quark-gluon shower and (d) conversion of the shower gluons into other 
quark antiquark pairs which (e) produce hadrons, followed by (f) the decay of the 
unstable hadrons. The hadronic showers carry the momentum of the original quark 
pair and so come out as two collimated "jets". Sometimes a hard gluon is also 
emitted in the process xx - qqg , which results in three jets. (More rarely two or 
more gluons can be produced.) The quark cascading and fragmentation processes must 
be described by quantum chromodynamics (QCD) models. The Lund simulation is based 
on a string fragmentation model of QCD interactions, with parameters' determined from 
the e e collider data. However, it enables one to take account of the different 
fragmentation effects associated with various quark species. It thus provides a 
framework for extrapolating from e + e~ collider data to determine xx annihilation 
spectra from differing mixtures of final state quark jets. 

Simulation of a xx annihilation begins by selecting a fermion-anti fermion final 
state according to the branching ratios (B.R.) appropriate for a given type of x 
particle which depend on the mass (m^) and charge (Qf) of the fermion (quark or 
lepton) involved (B.R. - BfFiif for photino (y) annihilation; B.R. - 8 f m f ^ for 
higgsino and Majorana neutrino (v M ) annihilations (Ref. 15, hereafter called" RS). 




a — - — — — c — d— - — e — f 


Fig. 1. Schematic diagram of the physics of x particle annihilation with the 
various stages of the secondary particle production described in the text. 
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ANTIPROTON PRODUCTION 


in e + e~ collider experiments, owing to 
However, e + e" annihilations produce a 
jets than that from xx annihilations 
states with a very rich mixture of b 
c quark jets carry off only 80% and 60% 

SRW and RS should be 
the xx annihilation 


The energy spectrum of cosmic-ray antiarotons from photino annihilation was 
first calculated by Stecker, Rudaz and Walsh (SRW), using p data obtained in e e 
collider experiments. They deduced from their results the ^fascinating possibility 
of explaining the entire observed spectrum of cosmic-ray p's as arising from the 
annihilation of 15 GeV x particles. RS extended the spectral calculations to 
antiprotons, x-rays and positrons, considering photinos, higgsinos and Majorana 
neutrinos as dark matter. Using the basic approach and formalism first given by SRW 
and RS, Stecker and Tylka 11 (ST1) extended the study of the anti proton spectrum with 
Lund Monte Carlo simulations, also taking account of recent upper limits on the low 
energy (< 0.5 GeV) cosmic-ray antiproton flux w ’ . 

SRW and RS argued that the spectra of antiprotons and y-rays would be quite 

similar to those observed for these products 
the universality of quark jet hadronization. 
different mixture of heavy and light quark 
(e.g., higgsino annihilation produces final 
quark jets). Collider data show that b and 
of the cms energy respectively. Thus, the results presented by 
considered to be a strict upper limit to the hardness of 

spectra. . . r 

As expected, the Lund simulation gives a significantly softer spectrum for 
final states weighted by heavy quark channels than that observed in e e 
experiments. Fig. 2 shows the inclusive antiproton production function in terms of 
the scaling variable x = E/M for higgsinos (and Majorana neutrinos) or_photinos of 
mass M obtained from the Lflnd simulation. E and 8 are the total p energy and 
velocity respectively. It should be noted that the cases of generic higgsinos and 
Majorana neutrinos of equal mass give identical results, as discussed by RS. Note 
that for x annihilations at rest the cms energy is 2M x> 

Fig. 2. Antiproton production 
spectra. Collider data are shown and 
compared with Lund prediction for 
e + e~ annihilation (solid curve). 
Lund predictions for 15 ^ GeV RFi 
(middle dashed curve) and yy (lower 
dashed curve) annihilation are also 
shown. The crosses show the Lund 
data points. 

The x annihilation spectra are 
compared with the spectral function 
obtained from the Lund simulation of 
e + e - annihilations at a collider cms 
energy of 35 GeV (ST1). The e e 
collider data are also shown. Fig. 
2 also shows the spectral source 
function used by RS to fit the TASSO 
data (upper dashed curve). The e e 
data do not include the antiprotons 
from n decay, however n decay was 
included in the xx annihilation 
calculations. For comparison with 
the e e and xx annihilation spectra 
in Fig. 2, the xx spectra were 
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multiplied by a factor of 0.5. ST1 fit the results from the Lund Monte Carlo runs 
to simple analytic functions of the form 

£pf( x )/B = (1/Bo) (do/dx) = A exp(-ax) + B(1-x) 6 /x (1) 

where c- is the number of p's per annihilation and f(x) normalized so that its 
integral is 1. This functional form reproduces the exponential dependence of the 
spectrum on x for low values of x. The second term dominates for x > 0.15 and goes 
smoothly to 0 at x = 1. The following source function fits to the parameter sets 
(A,a,B,s) in eq. (1) were found: (a) For the e + e _ collider simulation, (82, 46, 
0.21, 2.13), (b) for 15 GeV higgsinos or Majorana neutrinos, (148, 50, 0.34, 4.39), 
and for 15 GeV photinos, (101, 50, 0.20, 3.99). 

Photinos produce only -2/3 as many antiprotons as generic higgsinos or Majorana 
neutrinos (0.2 per annihilation as opposed to 0.3). Much of this difference is due 
to the t production channel, which is unimportant in higgsino annihilation but which 
accounts for 56% of photino annihilations and produces no baryons. (It should be 
noted that baryon-antibaryon production through quark jets is not as well understood 
as meson production (see discussion in ST1 and references therein.)) 

The production rate of antiprotons from xx annihilation in the halo is 

Q— ( E ) = n 2 <o v> f— (E) 12) 

P x x a p v ’ 

where f- (E) = dN-/dE, normalized to the number of antiprotons per annihilation. 
For example, fn the case of R's, the annihilation cross-section <av> is 
overwhelmingly dominated by the contributions of t leptons and c and b quarks in the 
final state and is 


<0 x V> A = 4~ ( + 3m c + 3m b ) = L3 x 10 " 26 cm3s_1 - (3) 


The antiprotons come from the hadronic cc and b5 final states, but these 
account for a fraction 6, = 30/31 => 1 of the total. The ratio of p yields from 

photino and generic higgstho annihilation is 


7 °p = K <0 y V> A /6 h <0 fi v> AJ 


= 9.4 x 10 ^(m,./m~) 4 <0.1 
' W q' 


( 4 ) 


where mu ^ 81 GeV is the W boson mass and m~ is the mass of the squark (scalar 
quark) that mediates the y annihilation (RS). 5 

The interstellar antiproton spectrum is of the form (from eq.(l)) 

1(E) = (4i)' 1 (p 2 /M 3 )<o v> CxB 2 {Aexp[-aE/M ] + BM [1-(E/M )] Y /E} (5) 

where p is the mass density of dark matter and t is the effective residence time of 
cosmic-fay antiprotons in the Galaxy. As pointed out by SRW and RS, eq. (5) for the 
interstellar spectrum must be modulated to take account of the effects of the solar 
wind on the spectrum observed at Earth. To do this, ST1 used a numerical 

integration of the spherically symmetric solution of the Fokker-Planck equation. 

Fig. 3 shows the resultant p/p ratio as a function of energy for higgsino and 
photino annihilation as compared with the present antiproton data, the calculated 
spectrum of antiprotons from cosmic-ray collisions and the possible spectrum of 
extragalactic primary antiprotons in the baryon symmetric cosmology, also 

modulated 11 (ST1) . For the higgsino case, the results are shown for both solar 
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Fiq 3 The ratio p/p as a function of kinetic energy given ror the different 
‘theoretical models (ST1). The data are as follows: (Ref. 17) I - black square; 
(Ref. 18) - black diamond; (Ref. 19) - black triangles; (Ref. 20). 


and solar 
cm. For 


minimum. The annihilation spectra 
generic higgsinos, <o.v> = 1.26 . 

a minimal squark mass of 80 teeV, giving 
from the figure that fluxes which are 


are shown for x 
x 10-^cmV 1 . 


<oV> = 0.22 
consistent 


maximum 
0.3 GeV 
case, ST took 
It is £l£ar 
limits 1B,iy on the 
anti protons in the 
These results 

time x is ~10 B y for halo propagation 
models . Studies of cosmic- 

indicate that x cannot be too much less than 2x10' y 
The density p is probably in the range ~0.2 to ~0.4 GeV cm 
higgsino value* for the cross section from eq. (3) is probably an upper 


= 10 8 y 

In t^s 

x 10" 26 
with 


p = 


and 

photi no 
cnrs" . 
the upper 


low energy antiprotons cannot account for the reported fluxes of 
energy range around 10 GeV. 
can be used to test for Restraints 


orux particles. The residence 

2xl0 7 y for (Lisle propagation 

•ray efectrons and secondary nuclei^’ appea n r 8 tn 
much less than 2xl0 7 y or too much 


and 


Rre 


to 

than 10 B y. 
The generic 
1 imit to what 


one would expect for a x annihilation cross section and significant values for the x 
masses lie in the range of -5 to -30 GeV (RS). For photinos the^cross sections can 


be much smaller than for higgsinos (see eq. (4)). Thus, r 5 = 10 


_ v- i 


(p/p) ?d is 
-4xl0^ 4 <ov> 


-10 (for 15 generic GeV higgsinos) but can easily be as low as 
The experimental upper limit r^ < 3.5 is consistent with an acceptable 
astrophysical parameters for all x particles of interest. I he 
supersymmetric parameters can be adjusted to give cosmological densities of 


at most 
<< 1 . 
range of 
relevant 
interest 


for values of M between -5 and -30 GeV. Generic higgsinos or Majorana neutrinos 
with M = 15 Getf will give o h‘ = 0.2. It is possible to generalize the values 

for th& vacuum expectation v&loes of the Higgs fields (see eq. (3) °f p S) 1 ° ne 
wishes to obtain a,h 2 = 1 (the inflation cosmology prediction) _ with 15 GeV 
higgsinos. Note theft 5 5ne requires a h 7 of at least -0.1 to obtain dark matter 
halos. (As usual, a is the fraction eff b Hie critical density in x particles and h 50 

is the Hubble constant in units of 50 km s Mpc _.) . 

The present upper limits on low energy p‘s do not constrain dark matter 
models. However, with the possible range of values for r 5 given above, it is 

possible to hunt for evidence of x's with low energy p experiments, perhaps at a 


level not too far below the present limits 
would place significant constraints on dark 


. Upper limits 
matter models. 


on 


the order of r^ < 0.1 
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THE COSMIC-RAY NEUTRINO SPECTRUM FROM X x ANNIHILATION 


There are four main sources of neutrinos and positrons from xx annihilation to 
consider, viz. , first generation prompt leptons, second generation prompt leptons 
and charged u and K meson decay leptons (RS). Charmed and bottom quarks and t 
leptons are efficient sources of prompt, high energy leptons and anti leptons from 
their weak decays. E.g. , for neutrino-positron production, the relevant decay 
chains are as follows (W* is a virtual W boson): 

"f* — . , 

T -+ V + W * 

T 

b” + c + W + * (6) 

and c + s + W + * 


with W * decaying to n v , where i = e or y (for b decays, a = x is considerably 
suppressed by phase space and will be neglected here). B.R.'s for these decays are 
18% for x and 13% for both b and c. Leptons also come from second generation prompt 
processes given by the decay chains x+y+e, c-*-y+e, b-^y+e and b-c+e. The B.R. for 
the first chain is 18%; that of the last three is 13%. 

The neutrino flux from the decay of tt mesons + produced in heavy fermion 
annihilations can be determined by using the e e” collider data. Within 
experimental error, the data for the pion production spectrum 
from e e annihilations with cms energies above 14 GeV can be fitted to a single 
spectral function having the form (RS). 


z f(E ) = b (11.6 e -1 ' 13E tt + 1.35 e~°- 46E 7r) GeV -1 . (7) 

7T TT 7T f V } 

2 - 21/2 

where, 6 = [l-(m /M ) x \ ' , z, is the pion multiplicity per annihilation. The 
low energy pion-3ec&y neutrinos, wRich peak at ~35 MeV, can be calculated using the 
well-known kinematical formulae. At energies » 35 MeV, the TT-decay neutrino 
spectrum may be approximated by noting that two pairs of v , x> and one pair of 
v , v are produced for each e , e - pair produced in th£ annihilations. All of 
tnese^leptons take about 1/4 each of the pion energy. It follows from eq. (7) that 


'(.) (E) 


<W e > Ct 

4 XT 


2 

n <o v> cx 

X X A 

4tt 


(5.4e- 1 - 84E ). 


(8) 


The muon neutrinos from kaon decay have a harder spectrum than the pion-decay 
component so that it gives a significant contribution to the neutrino spectrum at 
about 2 GeV. The kaon-decay neutrino spectrum as a function of v energy can be 
approximated by the expression y 




n <ov> cx 
x A 

4 TT 


(0.73e -0 * 76E ) 


(9) 


The energies of all light leptons produced in each of the other processes 
respectively are similar. Neutrino spectra of v 's (v 's) and v 's (v 's) may be 
calculated from the expression e e y y 

I(E v ) = Q(E v )<sl>/4tt (10) 

with <i> being the mean-path-length of the annihilation region ( e.g. through the 
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galactic halo). The integral spectra for M = 15 GeV are shown in Fig. 4 for a 
galactic center source. These spectra are normal ized to the upper limit on the y- 
ray flux from the galactic center assuming that y-rays at this level and. neutrinos 
could be from xx annihilation. i 

The differential v production spectrum has been calculated for 15 GeV 
higgsinos (or Majorana neutrinos) using a Lund model Monte Carlo program with the 
results given in Fig. 5. These results give a spectral shape which agrees well 
with the spectrum obtained from the analytic calculations given in Fig. 4. 



Fig. 4. The integral neutrino spectra from the galactic center normalized to the y- 
ray upper limit as discussed in the text. 

ANNIHILATION POSITRONS 

The positron source spectrum is similar to that for neutrinos. However, the cosmic- 
ray positron spectrum is changed in propagating to Earth ( RS ) because of of energy 
degredation of the positrons by Compton scattering and synchrotron radiation. A 

detailed discussion of this is given by Rudaz and Stecker (RS). The upper limits on 

the low energy p's imply that positrons from xx annihilation will have a flux which 
is expected to be significantly below that produced by cosmic ray interactions 
(ST1) . This- spectrum has been calculated using the Lund program by Tylka iJ and is 
shown in Fig. 6. This result is generally in good agreement with the spectral shape 
calculated semi-analytical ly by RS. However, the magnitude of the flux has been 
normalized taking account of the upper limits on low energy p's so that it is 

considerably smaller than that of RS, as well as the cosmic-ray induced flux. 

THE FLUX AND SPECTRUM OF NEUTRINOS FROM SOLAR CAPTURE AND ANNIHILATION 

28 

Galactic x particles can be captured by the gravitational field of the Sun 
and come together in the solar interior to annihilate, with the resulting flu^ of 
solar neutrinos potentially obervable by neutrino detectors on the Earth. T^g 
spectrum of solar annihilation neutrinos, which has been discussed by many authors 
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Fig. 5. The differential electron neutrino production spectrum from 15 GeV dark 
matter annihilation using the Lund Monte Carlo model 2 ' 
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Fig. 6. The cosmic ray positron spectrum from xx annihilation compared with that 
produced by cosmic ray interactions . Source spectrum - dashed Ti rip and right 
hand scale; spectrum with propagation over a mean lifetime of 3x10' y - solid 
line and left hand scale. 
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will be significantly different from the one shown in Fig. 4, owing to the fact that 
muons and tt's stop or interact in the solar interior before they decay. The 
resulting loss of v flux or energy eliminates the n and K meson decay components and 
most of the second generation components from the observable spectrum. Calculations 
with the Lund model give a softer spectrum than most of those g fq^nd previously, 
lowering the expected solar neutrino event rate by a factor of ~3. ’ 

Galactic x particles are captured and trapped by the Sun at a rate 

r tr - (3,/2) 1/2 (2G W , sc (11) 

where f is the probability that a x particle loses enough energy by elastic 
scattering to be captured ’ . Although lower mass x particles can evaporate from 

the Sun before annihilating, those of mass > 5 GeV, with which we are concerned, 
will not evaporate. Then, in equilibrium, the annihilation rate will be equal to 
half of the trapping rate and the neutrino flux at Earth will be given by 

I(E u ) = (r tr /2)(4.d 2 )- 1 t /(E v ) 

where d = 1 A. U. = 1.5 x lO^cm. 

For typical dark matter halo parameters n x 
the trapping rate is 

r tr - 10 65 [M x (GeV))- 1 Opx . px {cm Z ) s- 1 . 

For hiqqsinos and Mamrana neutrinos, the elastic scattering cross section in 
eq. (13), is o = 1.5 x 10~ 78 cm. Taking M x = 15 GeV, we find r tr - 10 s and 
the neutrino flux from solar xx annihilation will be 

l( x >(E ) = 2 x 10' 2 c f(E ). ( 14 ) 

v v V V 


10 1 


10 ° 


10 2 


10 3 


10 " 4 


10 5 


Fig. 7. The differential electron neutrino production function from 15 GeV X x 
annihilation in the solar interior. The curve marked LMC is from the Lund Monte 
Carlo calculations. 7 



E (GeV) 

v 


( 12 ) 

= 0.4/M x (GeV) and v = 300 km s"*, 

(13) 
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The total neutrino production function t f(E ) relevant to eq. (14) is shown 
in Fig. 7. Fig. 7 shows both the first g^ner v ation prompt (PI) approximation, 
similar ±o used in Refs. 32 and 33 and that obtained using the Lund Monte Carlo 
program , which gives a softer spectrum. Given the fluxes calculated using eq. 
(14), the event rate observed in a neutrino detector will be given by 


r ev - "n f dE v eff (E) K'UVU + • 05) 

where V gff is the effective detection volume, equal to the detector volume itself 
for the contained events from lower energy neutrinos (typically 0.5 to 2 GeV) and is 
proportional to E for through going detector events of higher energy . The 
neutrino-nucleon cross sections o and a- are linearly dependent on energy. Using 
eqs. (14) and (15), we can estimate the v ratio of solar annihilation to atmospheric 
events for the IMB neutrino detector, following the discussion of Ng, et al. . For 
contained events, this ratio is just 

2 2 

R C = ( 4 ^ 0 )r dE E I i x) < E >J / If dE E I A ™ (E) = 0.03 (16) 

L U L 0.5 v 0.5 v V ’ 

where is the solid angle around the Sun, determined by the detector resolution 
while for through going events it is 

R T = ( 4 "/«q)t fA dE E 2 ^ X) (E)] / [ j”dE E 2 I A ™ (Ej| = 1.5 (17) 

using, the production function from the Lund model. Fiq. 8 shows this production 
function along with the PI approximation and a typical E - ^ atmospheric spectrum (all 



Fig. 8. Relative event rate predictions as a function of energy for solar dark 
matter annihilation and typical atmospheric neutrino fluxes . Curve PI and LMC 
are as in Fig. 6. The atmospheric spectrum is typically an E“ 3 power law^. 
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weighted by E 2 as relevant to the throughgoing event rate,). This figure 
illustrates (1) that the event rate is overestimated by a factor of -3 by the simple 
first generation approximation, and (2) that the best "window" for observing solar 
annihilation neutrinos is the 2-10 GeV energy range. 

ANNIHILATION GAMMA RADIATION 

The spectrum of y-ray background radiation from X x annihilation in the halo may 
be calculated by noting that the continuum flux is overwhelmingly due to the decay 
of neutral pions produced in the xx annihilations. One can then ma^e use of the 
pion production spectrum (7) in order to determine the y-ray spectrum.. ^ 

The y-ray spectrum resulting from the decay of the neutral pions is given by 

M 

s f (E Y ) = 2 ; dE, (E 2 - mJ)- 1/2 c/(E „) (18) 

W 

where E (E ) = E + m 2 /4E and c is the y-ray multiplicity. 12 

A ifior£ exac£ calculation carS be made again using the Lund Monte Carlo program 
(ST2) . For x particles below the b quark threshhold (M = 5 GeV), -90 % of the 
resulting y-rays are from n decay. The remaining 10% come from the decay of other 
hadrons. For larger mass x particles, a component from B* meson decay produces a 
distinctively hard spectral signature in the -0.1 GeV energy range. This channel is 



Fig. 9. Gamma-ray spectral production functions for photino (open circles) and 
higgsino (closed circes) annihilation. The dashed longer curves show the pion 
decay component only. The 15 GeV plot shows the B* decay component 
from fifi annihilation separately. 
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especially important for R's which annihilate into a bB pair -85% of the time. In 
this case, y-ray production occurs -85% through pion decay, - 8 % through B* decay and 
-7% through other channels. Since y ' s have much smaller branching ratios for b 
production (-15%), B* decay accounts for only about 2% of the y-ray production in 
their annihilations. Fig. 9 shows the results obtained by Stecker and Tylka 12 (ST2) 
for y's (closed circles) and R's (or v 's)(open circles) with masses of 5, 6 and 15 
GeV expressed as the quantity 5 f(E ), where c is the average Y-ray multiplicity 
per annihilation and f(E ) is Y the Y normalized spectral distribution function. 
The n-decay component, is Y shown by a dashed line. The most striking feature is the 
contribution from B* -»■ B + y decay, particulary for Rfi annihilation. The decays 
B* -»■ B + y produce a y-ray line of energy £*- 51.7 MeV in the B* rest frame. For 
isotropic decays, this line transforms to a Square wave spectrum in the c.m. frame 
of the xx annihilation with sharp drops at the cutoff energies . For M = 15 GeV, 
the cutoff energies are 281.4 MeV and 9.5 MeV (see Fig 9). x 




Fig. 10. The high latitude galactic y-ray spectrum calculated for various 10 GeV x 
particle annihilations compared with the extragalactic (EG) diffuse background 
and galactic disk radiation (see text). 

Fig. 11 . The annihilation spectrum from a hypothetical dark matter core at the 
galactic center consisting of either 15 GeV y's or R's compared with the y-ray 
flux at the galactic center from cosmic ray interactions. 
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The high latitude galactic y-ray spectrum from xx annihilation is (RS) 


1(E) = 4^ > n 2 <ov> t f(E ) 
v y y 4 tt x &nn Y Y 

. 1.8 x 10‘ 6 <* 8 >»o.3 M x 2<oV> 26S f<E y ) 


-2 -1 -1 
cm s sr 


(19) 


where <2,> is the mean 1 ine-of-sight through the galactic halo in units of 8 kpc, 
Pn o is the x density in units of 0.3 GeV cnT J , and the annihilation cross section 
is in units of 10~ 76 cmV 1 . The fluxes at high galactic latitudes from the 

annihilation of 10 GeV R's and 10 GeV y's are shown in Fig. ? for pq 3 = 1 and <i> 
= 1. The generic R cross section is <oV> ?6 = 1.26. For y'Sj ST2 took a cross 
section corresponding to a lower limit on tne squark mass of = 80 GeV, which 
gives <ov> ?fi = 0.22. 

The annihilation fluxes are compared in Fig. 10 with an estimate of the mean 
diffuse background flux expected from cosmic-ray interactions with interstellar gas 
in the, gal art ic .disk at high galactic latitudes (mean total column density of 

-4x10^° cm -7 ) 17,38 . However, the high latitude interstellar medium is very 

patchy. Recent high resolution 21cm measurements, have shown that there are regions 

covering ~1% of the sky with neutral column densities < 0.9x10 cm (Ref. 39). In 

these regions, it may be possible to observe annihilation Y-radiat^o^from a naked 
dark matter halo using a high angular resolution Y-ray telescope. * 

There appears to be an extragalactic background of Y -radiation (ajry£ EG 
Fig. 10) which is softer than the galactic high latitude background * . The 

extragalactic flux may exceed the high latitude galactic background for energies 
below 0.4 GeV and may be comparable to the galactic component at higher energies 
provided that one can extrapolate from the observations. However, some theoretical 
models indicate that the extragalactic spectrum may steepen significantly above the 
observed energy range . 

In addition to a Y-Tay continuum, dark matter annihilation can produce line 
radiation in the GeV Y -ray region. For photinos of mass greater than 4 GeV and 
taking a minimal squark mass, Rudaz^ finds that the line flux in the direction of 
the galactic pole will be 


F — = 1 x 10 ^p 2 7 <Iq> (cm 2 s-sr) 1 (20) 

YY^YY U • o O 

independent of mass, with^the line centered around M and having a Doppler: width 
given by aE /E = 6 = 10 . It follows that if an erltergy resolution of 10” could 

be obtained! tKe 1 ifie-to-continuum ratio in a bin centered around the line would be 


F,. /F =< 0.7 (E /10 GeV) 1 * 7 . (21) 

line g v y 

Even with a 1% energy resolution, such lines should be detectable. 

It is possible for dark matter x particles to be concentrated in a core at the 
galactic center by the drag of ordinary baryonic matter through collapse in the 
early stages of galaxy formation^ . Annihilations from such a^source at a distance 
r $ consisting of x's in a volume V $ with mean-square density <p^> will give a flux 

F(E ) = ( 47 tr 2 ) _1 <p 2 >M" 2 <ov>V c f(E ) (cm 2 s) -1 . (22) 

v y s' xx s y Y 

Using the isothermal core model of Ipser and Sikivie^ 8 , we find 

(4nr 2 ) -1 V < P 2 > = 6.9 x 10 20 GeV 2 cnf 5 . (23) 

v s 7 s x 

For such a core of 15 GeV R's, the Y-ray flux is given by 
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F gc ( 0.15) - 1.5 x 10" 6 (cm^s-GeV) -1 (24) 

with a relatively flat spectrum between 0.1 and 0.2 GeV owing to B* decay. This 
spectrum is much harder than the y-ray spectrum from cosmic-ray interactions in the 
Galaxy, which acts as a foreground source. Both cosmic-ray (CR) induced and 15 GeV 
x annihilation spectra are shown in Fig. 11 for a y-ray telescope with a 1° beam 
size. The (CR) spectrum for the inner galaxy has been calculated theoretically for 
the inner galaxy and agrees well with the observational data . Let us consider the 
observablity of such a source with the EGRET y-ray telescope to be launched on the 
Gamma Ray Observatory (GRO).. This telescope has (1) an angular resolution of the 
order of 1° in the 0. 1-0.2 GeV energy range, (2) an effective area of -1500 cm 2 in 
this energy range, and (3) an energy resolution of -15% . Dividing up its 80-200 

MeV data into 3 energy bins of -40 MeV each, with a 1 month exposure time, even the 
lowest energy bin would produce a signal of -7 a. Using an on-off subtraction 

technique would further define the galactic center source. A detailed discussion of 
this work may be found ST2. 

The extragalactic and cosmological y-ray background spectrum from neutral heavy 
fermion annihilation can be calculated following the methods given by the author . 
and can be shown to be negligible compared to the observed extragalactic background. 

CONCLUSIONS 

We have discussed the annihilation physics of exotic dark matter particles, in 
particular, those predicted by supersymmetry theory. This physics leads to the 
production of calculable, and potentially observable fluxes of cosmic-ray and cosmic 
y-ray annihilation products. We conclude that a study of cosmic-ray antiprotons may 
give information about the nature of the dark matter, assuming that it is made up of 
x particles. In addition, studies of galactic y-radiation may also shed light on 
the dark matter problem. A characteristically hard spectrum in the 100 MeV region 
from a source at the galactic center or at high galactic latitudes could serve as a 
signature of x particle annihilation. With significantly more sensitive y-ray 
telescopes, the discovery of monochromatic radiation in the several GeV range would 
provide the most conclusive evidence. Solar neutrinos of several GeV energy may 
also give an observable signature for galactic x particle dark matter. However, the 
positron flux from dark matter annihilation would be buried below that from cosmic- 
ray produced secondaries. 
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